Acoustic Waveform Tomography: A Blind Test on A Synthetic Dataset

نویسنده

  • Chaiwoot Boonyasiriwat
چکیده

Acoustic waveform tomography is applied to a 2D synthetic dataset from Saudi Aramco to test if we will be able to obtain an accurate and highly resolved velocity tomogram. Source and receiver locations are not located along a flat surface so we assume a surface topography from the source and receiver locations. The initial velocity model was obtained using traveltime tomography and the source wavelet was estimated from the frequency spectra of the data. Two sets of lowfrequency data with peak frequencies of 2 Hz and 5 Hz were used in waveform inversion. Our waveform tomogram is not accurate compared to the true model. To investigate why waveform inversion failed, I generated synthetic traces using the true model, inverted them to get an accurate tomogram. Apparently some of the information about the model needs to be corrected so that the data file is reliable. After I obtain the correct acquisition information, I should get a more accurate inversion result. INTRODUCTION Waveform tomography inverts the full-waveform information in seismic data to get a high-resolution and highaccuracy velocity estimate of the subsurface geology. In practice, waveform tomography is a highly nonlinear inverse problem and very sensitive to the source wavelet, initial model and other factors. Some time-domain methods were proposed to reduce the nonlinearity of waveform tomography including early-arrival waveform tomography (Sheng et al., 2006) and multiscale waveform tomography (Bunks et al., 1995). Though waveform tomography has been developed for more than 20 years, it is still not a standard method for velocity estimation. Successful applications of waveform tomography are limited to a few examples. There are several factors that affects its success including an unknown source signature, choice of forward modeling (acoustic, elastic, or anisotropic wave equation), surface topography, correct selection of the number of model parameters, etc. Research is still needed to make waveform inversion a more robust and practical tool in exploration geophysics. In this work, I applied acoustic waveform inversion to a synthetic dataset generated by Saudi Aramco. The goal of this study is to test the sensitivity of our waveform inversion method to some factors mentioned earlier and to obtain an accurate velocity tomogram compared to the true model from Saudi Aramco. We have to deal with an unknown source wavelet, surface topography, and unknown forward modeling. The latter problem occurs when different forward modeling codes are used for synthetic trace generation and inversion. The waveform inversion results presented in this report were obtained using a rugged surface topography implied from the source and receiver locations. Later I discovered that the acquisition surface was flat. However, I do not have time to rerun the inversion using a flat surface. THEORY In this report, we consider wave propagation as described by the 2D acoustic wave equation, 1 κ(r) ∂p(r, t|rs) ∂t −∇ · [ 1 ρ(r) ∇p(r, t|rs) ] = s(r, t|rs), (1) where p(r, t|rs) is a pressure field at position r at time t from a source at rs; κ(r) and ρ(r) are the bulk modulus and density, respectively; and s(r, t|rs) is the source function. The forward solution, p(r, t; rs), of the twoway wave equation 1 is computed by a staggered-grid, explicit finite-difference method with 4-order accuracy in space and 2-order accuracy in time (Levander, 1988). The free-surface boundary condition is applied to the top boundary of the model, and the perfectly matched layer (PML) boundary conditions (Berenger, 1994; Chew and Liu, 1996; Zeng et al., 2001; Festa and Nielson, 2003) are 193 194 Boonyasiriwat utilized at the other boundaries. In terms of the Green’s function G(r, t|r, 0), one can write the solution p(r, t|rs) of equation 1 as

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Experimental Testing of the Blind Ocean Acoustic Tomography Concept

Acoustic focalization is a well known concept that aims at estimating source location through the adjustment of multiple environmental parameters. This paper uses the same concept for inverting water column sound speed in a blind fashion, where both source location and source emitted waveform are not known at the receiver that is Blind Ocean Acoustic Tomography (BOAT). The results obtained with...

متن کامل

A new stochastic 3D seismic inversion using direct sequential simulation and co-simulation in a genetic algorithm framework

Stochastic seismic inversion is a family of inversion algorithms in which the inverse solution was carried out using geostatistical simulation. In this work, a new 3D stochastic seismic inversion was developed in the MATLAB programming software. The proposed inversion algorithm is an iterative procedure that uses the principle of cross-over genetic algorithms as the global optimization techniqu...

متن کامل

Blind Ocean Acoustic Tomography: Experimental Results on the Intifante’00 Data Set

Blind Ocean Acoustic Tomography (BOAT) is an ocean remote exploration concept similar to acoustic tomography but where both the emitted signal waveform and the source position are unknown. BOAT consists of a minimal environmental model of the area, a broadband matched-field processor and a genetic algorithm search procedure. This paper presents the results obtained with BOAT on part of the data...

متن کامل

A Robust SAR NLFM Waveform Selection Based on the Total Quality Assessment Techniques

Design, simulation and optimal selection of cosine-linear frequency modulation waveform (CNLFM) based on correlated ambiguity function (AF) method for the purpose of Synthetic Aperture Radar (SAR) is done in this article. The selected optimum CNLFM waveform in contribution with other waveforms are applied directly into a SAR image formation algorithm (IFA) and their quality effects performance ...

متن کامل

Full-waveform inversion of GPR data acquired between boreholes in Rustrel carbonates

Full waveform inversion (FWI) of seismic or Ground Penetrating Radar data provides high-resolution quantitative images of the constitutive parameters of the rock/soil which control seismic/GPR wave propagation. We developed a 2D inversion tool in the frequency domain adapted to the multi-parameter physics controlling GPR propagation in isotropic non dispersive media, i.e. dielectric permittivit...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2009